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coefficients du predicteur soit en effectuant une transformee de Fourier sur les echantillons X(n). 




Xerox Copy Centre 



EP 0 410 870 A1 



DISPOSITIF D'ANALYSE SPECTRALE D'UN SIGNAL NUMERIQUE. 

La presents invention conceme un disposlttf d'analyse spectrale d'un signal numerique. 

Le brevet des Etats-Unis d'Amerique n * 4 787 055 decrit un tel dispositif qui utilise la methode connue 
sous le nom de Transformee de Fourier Discrete". Par cette methode on peut determrniner un spectre sur 
une large gamme de frequence. Cependant si Ton veut une analyse precise de ce spectre, on est amene a 
5 prevoir des temps de calculs qui croissent avec la precision demandee ; les temps de calculs peuvent alors 
devenir redhibltoires pour certaines applications oCi Ton cherche a obtenir une large gamme de frequence 
et une analyse fine du spectre. 

Une application visee est celle de Tanalyse d*un signal numerique elabore par un radioaitimetre a 
modulation lin§aire de frequence. On trouvera dans la demande de brevet frangais n* 2 600 778 la 
10 description d'un tel radioaitimetre. 

Dans ce genre d*application on est confronte a un premier probleme qui est celui de la grande variation 
d'aftitude a mesurer allant de quelques metres a quelques kilometres, et a un deuxieme probleme qui est 
celui de la rapidite de la mesure a fournir, compte tenu de la vitesse de Taeronef dans lequel est place le 
radioaitimetre- En outre les mesures doivent etre donnees avec une precision compatible avec la securite a 
75 assurer, surtout en phase d'atterrissage. . . » 

Le dispositif precite ne peut pas resoudre ce double probleme. 

La presente invention propose un dispositif de ce genre qui", lui, resout ce double probleme et est done 
bien adapte a Tapplication envisagee. 

Pour cela, un tel dispositif d'analyse spectrale est remarquable en ce qu'ii comporte des acces de 
20 ' commande pour recevoir une information de methode d'analyse destinee a commander au moins un circuit 
d'analyse. 

L'idee de rinvention repose sur le fait que Tanalyse spectrale doit etre effectuee selon la methode qui 
. est la mieux adaptee aux circonstances. 

Un dispositif d'analyse spectrale est remarquable, selon une caracteristique de rinvention, en ce que . 
25 Tune des methodes utilisees est la transformee de Fourier discrete, une autre impllqiie la methode du 
gradient et une troisieme la methode de Levinson. Ainsi. dans lecas de rappllcation envisagee a la mise en 
route du radioaitimetre la methode de transformee de Fourier discrete est mise en oeuvre pour analyser 
grossierement le spectre et si le spectre analyse laisse supposer une altitude elevee, la methode impliquant 
la rnethode de Levinson est utilisee et si I'altitude. est basse, c'est la methode impliquant la methode du 
30 gradient. ^ 

La descriptior) suivante,' accompagnee des dessins ci-annexes, le tout donne k titre d'exemple non 
limitatif, fera faien comprendre comment rinvention peut etre reaiisee. 
La figure 1 represente le schema d'un dispositif conforms a ('invention. 

La figure 2 montre le schema d'un premier circuit d'analyse faisant partie du dispositif montre a la figure 
35 1. 

La figure 3 montre la constitution des mots de commande du premier circuit d'analyse. 
La figure 4. est un diagramme-temps destine a I'explication de la formation des coefficients de 
correlation. • ' . 
La figure 5 represente I'implementation des memoires d'un circuit 6. 
40 La figure 6 est un diagramme-temps destine a rexplication de la methode de Levinson. 

La figure 7 represente I'implantation des memoires de circuit 6 lors de la mise en oeuvre de la methode 
de Levinson. 

La figure 8 est un diagramme-temps relatif a la methode de Levinson. 
La figure 9 est aussi un diagramme-temps relatif a la methode de Levinson. 
45 La figure 10 est un diagramme-temps destine a expliquer le fonctionnement pour la methode du 
gradient. \ ' 

La figure 1 1 represente I'implantation des memoires pour la methode de gradient. 

La figure 12 represente un diagramme-temps pour expiiquer le decienchement effectif de la methode du 

gradient. 

50 La figure 13 represente la structure du deuxieme circuit d'analyse. 

La figure 14 represente les mots de commande utilises pour ce deuxieme circuit d'analyse. 

La figure 15 est un diagramme-temps pour I'explication de la transformee de Fourier discrete. 

Le figure 16 est un diagramme-temps pour expliquer les operations efiectuees syr les coefficients d'un 

predicteur: 

La figure 1 represente le dispositif d'analyse spectrale conforme a rinvention. A la borne d'entree 1 il 
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regoit sous forme numerique les echantillons X(n) (oii n indlque I'instant d'apparition) d^un signal dent on se 
propose de faire Tanalyse spectrale. Ce signal est emmagasine dans un registre 2 sous la commande d'un ' 
signal EDB applique a la borne 3. Les composantes spectrales apparaissent sous fornie numerique, aussi. 
a la borne 4. 

6 Conformement a Pinvention le dispositil comporte des acces 5 pour recevoir une information de 
m^thode d'analyse destinee k commander les circuits d'analyse 6 et 7. 

Dans rexemple decrit, trois methodes sont envisagees. mals ceci ne limite en rien I'invention. 
Ces methodes envisagees sont : 
La methode de la transformee de Fourier discrete executee uniquement par le circuit 7, 
10 Les methodes du gradient et de Levinson qui, executees par le circuit 6. fournissent les coefficients d'un 
predicteur au circuit d'analyse 7. Ce circuit 7 etablit finalement les composantes spectrales. 

Les methodes de Levinson et du gradient fournissent les coefficients a(k) d'un polynome predicteur qui 
determine un echantillon predit XP(n) en fonction des echantillons reels regus precedemment ; 

XP(n) = E a(p)X(n-p) (1) 
P=1 

ou M definit Tordre de modelisation. 
20 A parttr des coefficients a(p) on etablit les composantes spectrales S(k) : 
S(k) = 1/{[DR(k)F + [Dl(k)p}* (2) 

DR(k) = 1 + E a(i).cos2irk.J (3) 

DI(k) = E a(i).sin2irk i (4) 
i=1 " 

1/N definit I'espace qui separe les composantes spectrales ; pour definir cet espace on prendra plutot PR ' 
= 1/N ; 
k = 1....NBR 

ou NBR definit le nombre de raies possibles et k leur rang. 

La methode de transformee de Fourier discrete fournit aussi les composantes S(k) : 
S(k) = {[SR(k)F + [SI(k)F}- (5) 

SR(k) = 1 + E X(i).cos2irk,^ (6) . 

i=1 

SI(k) = E X(i).,sin2irk ^ (7) 

On notera aussi que !a methode de Levinson fait intervenir les coefficients de correlajion de la suite 
d'echantillons X(n) ; 

1 NN-1 \ny 

R(j) = 15]J E X(n).X(n+j) (8) 



30 
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50 

OLf j = 0 ... M + 1 

et NN represente la longueur, ou horizon d'accumuiation. On fait aussi intervenir un parametre d'oubli L 
dans la methode du gradient pour un ajustement du coefficient d'autocorrelation R(0) a chaque echantillon 

X(n): 

55 R(0) = (1-L).R(o) + X(n)2 (9) 
onaO^(l-L)^1 

Les parametres definissant ces operations, notamment les vaieurs M, PR, NBR. NH, L sent transmis a 
racces 5 de la maniere suivante. 
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L'acces 5 se decompose en une ligne de donnees 10 et une ligne d'adressage 13. La ligne 
d*adressage reliee a un circuit de regi'stres de parametres 15 permet de selectionner un des registres 21, 
22. 23. 24, 25 et 26 contenu dans ce circuit 15 de sorte que les donnees presentes sur la ligne 10 / sofent 
enregistrees. 

5 Les registres 21 et .22 sent destines a contenir les informations regissant Texecution de la methode du 
gradient et de Levinson c'est-§t-dire. entre autres, le parametre d'oubli L et I'ordre de modelisation M. 

Les registres 23 a 26 contiennent des informations concemant le fonctionnement du circuit d'analyse 7, 
notamment les parametres PR, NBR... 

Ceci sera expiique plus en detail dans la suite du present memdre. 
10 Le dispositif d'analyse comporte egalement une horloge de datation 30 servant a foumir une date a des 
evenements, horioge couplee a un registre de datation 31 pour emmagasiner des dates remarquables. 

Un circuit de sortie 35 penmet d'aiguiller sur la sortie 4, outre les composantes de spectre, la valeur de 
coefficient d'autocorrelation R(0) souvent utile pour reguler le niveau d'amplification.d'une chatne reception, 
le contenu du registre 31 et la date elaboree par I'horioge 30, plus une valeur d'un statut definissant I'Stat 
IS dans lequel se trouve le dispositif- Differents amplificateurs a trois etats 38 et 39 permettent a une seule 
ligne commune 40 de transmettre, sur races 4f les informations de sortie du circuit 35 et les infonmations 
d'entree pour le circuit 15. L'amplificateur 39 est mis a I'etat passant par un signal RWB exterieur au circuit 
et mis a I'etat actif pour signaler qu'une information doit etre enregistree dans (e circuit 15. L'amplificateur 
38 est mis a Tetat passant par un signal TS emanant du circuit 15 en reponse a des demandes de 
20 Texterieur appliquees a racces 5. Une horloge de cadencement 42 pilotee a partir de I'ext^rieur fournit des 
signaux H pour fixer la cadence de fonctionnement de tout le dispositif. 
Differents signaux sont appliques au dispositif de Tinvention : 

- le signallNI pour initialiser tout ie dispositif, 

- le signal RZ1 pour initialiser le circuit 6. 
25 - le signal R22 pour initialiser le circuit 7, 

- le signal INHB pour permettre certaines operations comme cela sera indique plus loin. 

On decrit dans ce qui suit, plus en detail, la structure et le fonctionnement du dispositif d'analyse 
spectrale de I'invention. 

1 - DESCRIPTION DU CIRCUIT D'ANALYSE 6 (BLOC 1) 

II se compose essentiellenrient d'un circuit sequentiel 50 et d'un circuit de calcul 51. 

Ce circuit sequentiel fournit differents signaux pour que le circuit de calcul 51 execute les differentes 

35 taches qui lui sont affectees. . . 

Le circuit de calcul 51 traite les echahtiilons numeriques qui sont appliques a sa borne d'entree 55 et il 
fournit a ses sorties 57 et 58 les coefficients a{p) et le coefficient de correlation R(0). L* analyse des mots 
MOT1 et M0T2 contenus respectivement dans les registres 21 et 22 permet au circuit sequentiel 50 de 
piloter le circuit de calcul pour que la methode de Levinson ou du gradient soit executee. Le circuit de 

40 . calcul est forme d*un ensemble de memoires BI, 62 et 63, d'un organe de multiplication 65 et d*une unite 
arithmetique et logique 67; d'un ensemble de registres 70, 71. 72, 73, 74 et 75 qui autorisent un 
fonctionnement en "pipe-line", d*un ensemble de multiplexeurs 79, 80, 81, 82 et 83 pour permettre 
I'aiguillage des informations vers les differents elements constitutifs du circuit 51, d'un organe de decalage 
85 et d*un circuit d'aide a la division 87 cooperant avec Tunite 67. Ces differents elements sont connectes 

45 de la maniere suivante. La sortie du multiplexeur 79 muni de deux entrees est reliee a Tentree de donnees 
. de la memoire 61 Tune de ses entrees constitue I'entree 55, Tautre entree est reliee a la sortie du registre 
75 qui constitue la sortie 57. La sortie du registre 75 est. aussi. reliee aux entrees des memoires 62 et 63. 
Les registres 70 et 71 sont destines a contenir les operandes pour la multiplication a effectuer par Torgane 
de multiplication 65, Les multiplexeurs 80 et 81 a quatre entrees correspondant a des positions de 

50 multiplexage permettent le stockage de ces operandes. Dans une position de multiplexage du multiplexeur 
80, ie registre 70 est boucle sur lui-m§me, dans une deuxleme position le registre 70 est mis en relation 
avec le circuit 87, dans une troisieme position il est mis en relation avec la sortie de donnees de la 
memoire 61 et dans une quatrieme position avec la sortie de donnees de la memoire 63. Quant aux 
positions de muftiplexage du multiplexeur 81 elles mettent en relation le registre 71 avec la sortie de 

55 donnees de la memoire 61 pour une premiere position, avec la sortie de donnees de la memoire 62 pour 
une deu>:rsme position, avec la sortie de Tunite 67, via rorgane de decalage 85 pour une troisieme position 
et avec lui-meme pour ia quatrieme position. Le multiplexeur 82 est a cinq positions et il met en relation le 
registre 72, qui est un registre d'operande pour i'unite 67, avec lui-meme pour une premiere position, avec 
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la sortie du registre 74 connectee eri sortie de Tunite 67 pour une deuxieme position, avec la sortie de 
rorgane de multiplication 65 pour une troisieme, avec la sortie de donnees de la mimoire 63 pour une 
quatrieme et avec la sortie de rorgane 85 pour la demiere ou cinquieme. Le multipiexeur 83 est a quatre 
positions et met en relation le registre 73, qui est un registre pour i'autoe operands de I'unite 67 avec la 

5 sortie d© donnees de la memoire 61 pour une premiere position, avec la sortie de donnSes de la memoire 
63 pour une deuxifeme position, avec le registre 76 dont I'entree est connectee a la sortie de I'organe de 
decalage 85 pour une troisifeme position et avec cette sortie de I'organe 85 pour une quatrieme position. 

Le circuit s^quentiel 50 foumit diffirents signaux SX1 . SMA. SMB et SAA pour I'aiguHlage & provoquer 
a I'interieur des multiplexeurs 79, 80. 81, 82 et 83 respectivement. Des codes de commande OPAL pour 

10 runU 67 peuvent §tre accompagnes de codes de decalage SHAL agissant sur I'operande contenu dans le 
registre 72. Ce drcuft s6quentiel 50 foumit aussi des codes d'adressages de lecture et d'6criture pour ies 
memoires 61, 82 et 63. Soient Ies codes ADXIW et ADX1R pour I'adressage en ecriture et en lecture de la 
memoire 61. Ies codes ADX2W et ADX2R concernant la memoire 62 et Ies codes ADRW et ADRR pour !a 
mgmoire 63. Ces codes d'adressage d'&riture et de lecture sent appliques a des multiplexeurs d'adressa- 

15 ge 91. 92 et 93 qui foumissent Ies codes d'adressages ADX1. ADX2 et ADR pour Ies memoires 61 a 63. 
Les commandos de commutation a I'interieur des multiplexeurs 91 a 93 sont effectu6es syst^matiquement 
au rythme du signal H. Une premiere demi-periode de ce signal est affectee a I'ecriture done, les codes 
ADX1 ADX2 et ADRR sont les codes ADXIW, ADX2W et ADRW. Pendant cette premiere demi-periode, 
une inscription pourra avoir lieu si le signal d'ecriture WX1. WX2 et WR applique aux memoires 61, 62 et 

20 63 est actif . 

2 - FONCTIONNEMPJT DU CIRCUIT D'ANALYSE 6 

25 • II est declenche par le decodage de I'ensemble des mots MOTI et M0T2 contenu dans les registres 
21 et 22 a huit elements binaires. La structure de ces mots est representee, a la figure 3. Un element 
binaire AUTOCO : la valeur "I" de cet element binaire declenche le calcul des coefficients d'autocorrelatfon 
et implique que la methode de Levinson va etre executee. Des elements binaires definissent M, qui fixe 
I'ordre de modelisation, I'element L fixe la valeur d'oubli utilisee dans ia methode du gradient ; AUTOB. acttf 

30 lorsque Pelement binaire est egal a "0", commande une remise a z6ro automatique et DRC considere pour 
la methode du gradient donne le depart de r^actualisation des coefficients, actif pour une valeur "1" de 
I'element binaire. 

35 2-A METHODE DE LEVINSON - A1 - Calcul des coefficients d'autocorrelation 

On se reporte a la figure 4 qui est un'diagramme-temps et a la figure 5 qui montre I'implantation dans 
les memoires 61 et 63 des differents valeurs en fonction des codes d'adresse. Pour la simplicity des 
explications on a limits a 5 le nombre de code d'adresse mais la generalisation a "p" codes se fera 
40 facilement. 

Sur les diagrammes-temps mentionnes dans le present expose, on a fait figurer des cycles de 
fonctionnement dont la duree est egale a la periode du signal d'horloge H. Conformement a ce qui a ete 
dit. cheque cycle est divise en deux parties. Durant la premiere partie. des ecritures peuvent avoir lieu dans 
les differentes memoires et la deuxifeme partie est consacree a la lecture. 

45 A ia figure 4, on se place a I'instant tO ou un echantillon valide surgit a la borne 55, soit X5 cet 

echantillon. La figure 5 montre I'implantation des echantillons antsrieurs X4, X3 XO dans la memoire 61 

ranges respectivement aux adreses •'0^ "1" .... "4". A cet instant le code d'adresse de lecture et d'ecnture 
de la memoire 61 est "4". I'echantillon X5 est inscrit dans cette memoire puis est lu et enfin se trouve 
enregistre dans les deux registres d'operande 70 et 71 de I'organe de multiplication 65 a I'instant t1 de 

50 sorte que le produit X5.X5 se retrouve dans le registre 72 a I'instant 12, tout comme le coefficient R(0) se 
retiouve dans le registre 73, ce code d'adresse de lecture de la memoire a l instant anterieur etant "0". La 
somme X5.X5 + R(0> se retrouve dans le registre 75 a I'instant t3 et concerne le code d'ecriture de la 
memoire 62, "0" ; le coefficient de correlation R(0) est done remis a jour par la contribution de X5.X5. La 
mise a jour du coefficient R(1) commence a I'instant t1 ou I'echantillon X4 est adresse en lecture et se 

55 retrouve a I'instant t2 dans le registre 71 et, comme le registre 70 contient toujours I'echantillon X5, le 
produit X5.X4 est forme a I'instant t3 dans te registre 72. Le coefficient R(1) se retrouve dans le registre 73 
pour y etre additionne avec le contenu du registre 72. le coefficient R{1) est re-enregistre dans la memoire 
63 a I'instant t4. Ainsi, de la m§me manifere, les differents coefficients restarts R(2). R(3) et R(4) vont etre 
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remis a jour par la contribution de I'echantiHon X5. Le processus se temnine a IMnstant t7. A Tinstant t8 
surgit un nouvel echantillon X6. Le traitement s'effectue alors de la meme maniere que pour Techantillon 
X5. On notera cependant le decalage des codes d'adresses d'ecriture et de lecture de la memoire 61. 

Le processus d 'autocorrelation s'arrete dhs que rhorizon programme est atteint La methode de 
Levinson est alors executee. 



2-A MCTHODE DE LEVINSON - A2 - Phase d'iniUaiisation 

10 Durant une premiere partie de cette phase, les coefficients du predicteur a(k). honmis a(0) = 1, sont 
mis a zero. La figure 6 est un diagramme-temps qui reprSsente cette phase d'initialisation ^ILC. Les codes 
d'ordre d'adressage ADF0<1(R) et ADRX2(R) vont evoiuer entre les valours "O" et "M-l", la valeur M etant 
contenue dans le registre 21. . * . 

L'evofutibn de ces codes n'a pas a etre consideree pour cette phase. On considere maintenant j'instant 

15 til oD le code OPAL applique a Tunite 67 est "0", c'est-a-dire que le resultat a sa sortie est a la valeur "0" 
quelles que soient les valeurs des operandes appliquees a ses entrees. Cette valeur "0" est enregistree 
dans le registre 75 et est enregistree dans les memolres 61 et 62 a Tadresse "4" apres Tinstant t12. A cet 
instant t12, le code OPAL est "1". ce qui impiique que le code a sa sortie est forcement "1", pour qu'a 
. I'instant t13, le chiffre "1" sort enregistre aux adresses "3" et "0" des memoires 61 et 62. Apres cet instant 

20 t13, le code OPAL redevient egal a "0" pour que des "0" soient enregistres dans les autres emplacements 
des memoires 61 et 62, Geci se termine a I'instant t15. La repartition des contenus des memoires 61. 62 et 
63 est montree a la figure 7. 

Durant la seconde partie de cette phase d* initialisation un parametre K intervenant dans la methode de 
Levinson est initialise, pendant une phase ^DIV a une valeur telle que : 

25 K = R(1)/R(0) en general K = V/W 

Cette division fait intervenir Torgane d'alde a la division 87 qui est un registre a decalage qui contient le 
resultat de cette division en emmagasinant au fur et a mesure les elements binaires constituant ce resultat. 
Cette division commence a Tinstant t15. Le dividende est contenu dans le registre 72 dans lequel a deja 
ete enregistree la valeur. R(1) a I'instant t13 et le divlseur est contenu dans le registre 73 qui contient la 

30 valeur R(0) depuis Tinstant t12. La comparaison des contenus des registres 72 et 73 donne la valeur d'un 
element binaire du quotient, et le contenu du registre 72 s'impute, en valeur absolue, de la valeur du 
diviseur qui reste identique a lui-meme durant tout le processus de division mais le contenu du registre 72 
est multiplie par 2. A I'instant t15 le code OPAL est -A.+ B ou -A-B selon le signe du dividende et du 
diviseur de du dividende et du diviseur de sorte que la quantite -K est en fait calculee. 

35 Cette valeur se trouve a Tinstant t16 apres le processus de division qui a dure 16 coups d'horloge 
correspondant a 16 bits, par exemple. . . 
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2-A METHODE DE LEVINSON - A3 - Phase de reactualisation des coefficients 



Cette phase impiique M fois I'execution successive de trois phases elementaires la phase : ^DIV deja 
mentionnee. la phase <^ACL qui modifie les coefficients, la phase <#>WL inseree dans la phase <#>ACL et 
prevue pour modifier un parametre W et la phase <>KL qui modifie la valeur des coefficients pour la phase 
<>ACL suivante. La valeur a(0) reste egale a "1" durant tout le deroulement de Texecution de la methode de 
45 Levinson. , 

La figure 8 est un diagramme-temps qui represente plus particulierem'ent la phase <»ACL. Apres la 
phase d'inltialisation. les operations suivantes sont executees durant la phase <^ACL. 
a(1)2 = a(1)i + K.a{0), = K (A1) 
a(2)2 = a{2)i + K.a(4)i = 0 (A2) 
50 . a(3)2 = a(3)i + K.a(3)i = C (A3) 

a(4)2 = a(4)i + K.a(2)i = 0 (A4) ou les indices 1 et 2 representent respectivement les memoires 61 
et 62. La premiere fols la memoire 61 est la source des coefficients a(i) et la merrloire 62 la destination 
de ces memes coefficients. 

A I'instant t16 (voir figure 8), le code ADRX2(R) est "0", ce qui permet au coefficient a(0) de se 
55 retrouver a I'instant tl7 dans le registre 71 et comme le parametre -K est dans le registre 70, le produit PC 
= K.a(O) se trouve, a I'instant t18, dans le registre 72 alors que dans le registre 73 se trouve le coefficient 
a(l)i et comme le code OPAL est -A + 8 le resultat a(1) + K.a(0) est tourni par l!unite 67 pour etre range a 
I'instant t19 eri transitant par le registre 75 dans la memoire 62. le code ADRX2(W) etant "1". Ainsi 
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roperation (A1) est executee. A cet instant t19 le coefficient a(4) provenant toujours de la memoire 62 se 
trouve dans le registre 71 pour que le produit P(4) = K.a(4) soit effectue. Ce produit est additionne avec le 
coefficient a(2) a Tinstant t20 de sorte que Toperation (A2) est executee et le nouveau coefficient a(2) est 
range dans . la memoire 61 a Tinstant t21. D'une maniere semblable les operations (A3). (A4) sont executees 
5 pour que les coefficients a(3) et a(4) soient ranges dans la memoire 62 aux instants t23 et t25. 

Du fait du fonctionnement en "pipe-line" on notera la presence d'une phase <>WL (voir figure 8) qui 
demarre avant la iin dela phase <J>D1V. Cette phase se deroule entre les instants t16 et t17 et a pour but de 
remettre a jour un coefficient w' a partir de sa valeur anterieure et du parametre K selon la fomnule : 
W' = W - K.V 

70 La valeur de ce parametre W est contenue dans un registre 76. Sa valeur initial© est : 
W = R(0) V 

Comme a Tinstant t16 on a obtenu la valeur -K et que la valeur V, dont la valeur initiale est R(1). est 
stockee dans le registre 71 , ie produit K.V est execute et se retrouve a Tinstant t17 dans le registre 72 pour 
etre additionne a la valeur W precedente provenant. du registre 76. La. nouvelle valeur w' est finalement 
75 emmagasinee dans le registre 76 a Tinstant t18. 

A la suite de la phase <;>ACL se deroule une phase ^KL dont le diagramme-temps est montre a la figure 

9. 

Durant cette phase, Toperation suivante est executee : . 
. V = R(2).a2(0) + R(1 ).a2(1 ) + R(0).a2(2)-R(5).a2(3) + R(4).a2(4) pratiquement : 
20 V = R(2) + R(1).a2(1) 

.:puisque les autres coefficients aaC) sont encore nuls. 

Cette phase debute a Tinstant t24 lorsque le code ADRX2(R) est le code "0" et comnrie le code ADRR- 
(R) etait deja egal a 2, les valeurs R(2) et 32(0) peuvent etre mises a IMnstant t25 dans les registres 70 et 
71. Le produit PR(0) de ces valeurs est a Tinstant t26 dans le registre 72. puis a Tinstant t27 dans le registre 
25 73 en transltant par Tunite 67. A I'instant t27 le produit PR(1) qui est le produit de R(1) par 82(1) se trouve 
dans le registre 72 de sorte que la somme peut §tre effectuee. le resultat SOI est alors mis, a Tinstant t29, 
dans le registre 73 qui va faire ainsi office d'accumulateur. Les produits successifs PR(2). PR(3) et PR(4) 
egaux respectivement a R(0).a2(2). R{5},a2(3), et R(4).a2(4) vont se trouver successivement aux instants t28, 
t29 et t30 dans le registre 72 et s'accumuler dans le registre 73 et le resultat final stocke a t31 dans le 
30 registre 71 et dans le registre 72. / , 

Apres IMnstant t31, une autre phase de division est declenchee pour evaluer la quantite V W ce qui 
donne le coefficient K utile pour revolution des coefficients du predicteur. Avant d'exarniner la phase 6ACL 
qui suit et qui debute a I'instant t32 .il convient de noter la valeur des codes d'adresse de lecture des 
memoires 61 . 62 et 63 (le code d'adresse d'ecriture etant retarde des codes d'adresses de lecture de deux 
35 temps d'horloge TH) qui sont respectivement (1 , 0, 3). 

A la phase suivante <j>ACL les operations effectuees sont ; 
a(1)i = a(1)2 + K.a(1)2 = K + K K (AS) 
a(2), = a(2)2 + K'.a(0)2 = k' (A6) 
a(3)i = a(3)2 + K'.a(4)2 = 0 (A7) 
40 a(4)i = a(4)2 + K ,a(3)2 = 0 . (A8) . 

Pour ces operations le role des memoires 61 et 62 est echange. de sorte que la memoire 62 va etre 
adressee comme la memoire 61 et vice versa, 

Ainsi, apres I'execution de plusieurs phases successives, revolution de la valeur des coefficients du 
predicteur est la suivante. donnee par le jeu d*equations (A9 a A12 d'une part et A13 a A16 d'autre part) : 



a(1)2 




a(l)i 


+ 


K .a(2) 


(A9) 


a(2)2 




a(2), 




K"a(1) 


(A10), 


a(3)2 




a{3), 


+ 


K".a(0) = 


k" (All) 


a(4)2 




a(4), 


+ 


K>(4) = 


0 (A12) 


a(1)i 




a(1)2 


+ 


K .a(3) 


(A13) 


a(2), 




a(2)2 


+ 


K'.a(2) 


(A14) 


a(3), 




a(3)2 


-f 


K".a(1 } 


(A15) 


a(4), 




a(4)2 


+ 


X".a(0) = 


k" {A16) 



55 2-B METHODE DU GRADIENT 

Cette methode va etre decrite pour un ordre 5, toujours a titre d'exemple. 

Elle est executee en quatre phases : une phase d'attente d'echantillon 6AG, une phase de calcul de R- 
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(0) : (#>RG, une phase de calcul de Terreur it>EG et une phase de remise a jour des coefficients du 
predicteur 4^CG: Cette methode est explicitee a Taide du diagramme-temps de la figure 10 et de 
rimplementation des memoires representee a la figure 11. 

5 * • • 

2-B METHODE DU GRADIENT - B1 - Phase 4f>AG 

Soit nnstant t50 ou apparaH un echantillon X(n). La phase d'attente se situe done anterieurement 
a cet Instant. Cette phase se caracterise essentiellement par un arret de revolution des codes d'adressages 
10 ADF0(1 et ADRX2. 



2-B METHODE DU GRADIENT - 82 - Reactuallsation de R(0) (phase *RG) 

75 Durant cette phase, la valeur de R(0) qui est contenue. dans le:registre 74 est remise a jour selon la 
formule (3) : 
R(0) = (1-L)R0 + X2(n) 

. La valeur L provient du registre 21 . Pour signaler qu'un echantillon X(n) se trouve a Tentree de la memoire 
61 via le multiplexeur 79, une transition descendante d'un signal EDI est . engendree. La valeur R(0) 
20 anterieure, contenue dans le registre 74 est mise dans le registre 72. Le code ^operation applique a I'unite 
67 est un code de decalage d*operande"As4" pour diviser Toperande contenu dans le registre A d'un 
multiple de deux, ce qui revient a le decaler d'une quantite depiendant de "L". Ce code decale se retrouve 
a'rinstant t51 dans le registre 73. Dans le registre 72. se trouve toujours la valeur anterieure R(0). 

A iinstant t50 rechantiilon X(n) est enregistre dans la memoire 61 a I'adresse "5" et prend ainsi la 
25 place de rechantiilon le plus ancien (a la figure 11, rechantiilon X(n-6)) et est restitue a I'instant t51 pour 
etre emmagasine dans les registres 70 et 71. Le produit X2(n) est alors execute. 

A rinstant t52, d'une part la valeur R(0) a subi I'influence du facteur d'oubli et porte la reference R°(0) 
-a la figure 10 cette valeur est emmagasinee dans le registre 74- d'autre part le resultat du produit X2(n).est 
mis dans.le registre 72. La phase ^RG est terminee a Tinstant t53, lorsque Taddition des valeurs contenues 
30 dans les registres 72 et 73 est effectuee et le resultat range dans le registre 74. 

2-B METHODE DU GRADIENT • B3 - Calcul de I'erreur (Phase <»EG) 

35 Durant cette phase le calcul suivant va etre execute : 

■ ■ -• 5 . . • 

e{n) = E ao.X(n-p) (B1) 

40 

ou e(n) est Terreur. 

. Cette phase demarre avec les calculs deja commences durant la phase ^RG. Ainsi a rinstant t5l les 
codes d*adresses de lectures ADRX1 (R) et AbRX2(R) etaient "5" et "0" respectivement, ce qui permetaux 
valeurs X(n) et a(0) de se retrouver dans les registres d^entree 70 et 71 de Torgane de multiplication a 
45 rinstant t52 de sorte qu'a I'instant t53 le premier element de la somme indiquee dans la derniere formule 
. se trouve dans le registre 72 a Tinstant t53 et a I'instant t54 ce resultat se retrouve dans le registre 73 qui 
va servir d'accumulateur. Ainsi aux instants t54. t55, t56, t57 et t58, les differents produits al.X(n-l ), a2.X(n- 
2). a3.X(n-3), a4.X(n-4) et a5.X(n-5) vont se trouver successivement dans le registre 72 pour que 
raccumulation puisse s'effectuer dans le registre 73. A rinstant t59. le calcul de e(n) est termine et le 
50 resultat est mis dans le registre 70 pour que la phase suivante puisse s'effectuer. 

2-B METHODE DU GRADIENT - B4 - Mise a jour des coefficients phase <)CG 

S5 Durant cette phase les calculs suivants sont effectues : . - . 

a'(0) = 1 (E2) 

a (p) = a{p) - TTT.e(n).X(n-p) (B3) 
avec p= 1 a 5, • 
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les valeurs accentuees representent la valeur mise a jour, la valeur TTT est normalisee par la valeur R(0) 
lors d'un calcul prealable. est fixe et est obtenue par un decalage de I'operande contenu dans le registre 
72. 

Ce code represente a la figure 10 par "-Ash'+ B" apparatt des I'instant t60, pour le calcul des differents 
5 coefficients de la formule (B3). La formule (82) etant calculee a Tinstant 159 le code d'operation applique a 
I'organe force a "1 " sa valeur de sortie. A I'instant t66 tous les coefficients contenus dans la memoire 62 
ont 4te mis a jour. 

10 2-B METHOOE DU GRADIENT - B5 - EnchaTnement avec un nouvel echantiilon 

A la phase <^CG peut se superpbser une phase d'attente pour Techantillon suivant X(n + 1). Sur la figure 
on admet que cet echantiilon surgit a I'instant t66. Get echantiilon est enregistre a Tadresse "4" de la 
memoire, le traitement s'effectue connme ci-dessus mais tous les codes d'adresses sont incrementes d'une 
IS unite et ainsi de suite pour les autres echantillons qui surviennent. 

II va de soi que differentes mesures ont ete prises pour assurer la coherence du calcul de coefficients, 
lorsque le dispositif est utilise dans un radioaltimetre du type d4crit dans la demande de. brevet 2622021. 

20 3 - DESCRIPTION DU CIRCUIT P'ANALYSE 7 (BLOC 2) 

' Le schema est represente a la figure 13. Ce circuit est constrtue essentiellennent par un circuit 
sequentiel 100 et par un circuit de calcul 101. 

Tout comme pour le circuit 5. le circuit sequentiel 100 fournit differents signaux de commande pour que 
25 le circuit de calcul 101 execute les. differentes taches qui lui sont affectees en fonction des informations 
contenues dans les registres 23. 24. 25 et 26 eux-m§mes contenus dans le circuit de registres de 
parametres 15 (figure 1). 

Le circuit d'analyse est muni de deux entrees 102 et 103 pour pouvoir traiter soit les echantillons pour 
effectuer une transformee de Fourier discrete, soit les. coefficients du predicteur afin de restituer les 

30 composantes spectrales. 

Les echantillons et les coefficients sont enregistres dans une memoire 105 en transitant par un 
multiplexeur 107. II y est prevu un organe de multiplication 110 pour faire le produit des nombres contenus 
dans des registres de "pipeline" 112 et 113. Ces registres sont connectes en sortie de la memoire 105 et 
d'une memoire 115 fournissant des valeurs de cosinus et de sinus, comme cela sera explicite dans la suite 

35 du present expose. Une unite arithmetique et logique 120 effectue differentes operations en fonction d'un 
code d'operation 0PA2. Les registres d'operandes 122 et 123 sont connectes respectivement aux sorties de 
multiplexeurs 124 et 125. Le multiplexeur 124 est muni de deux entrees ; I'une regoit un code d'initialisation 
pour initialiser le contenu du registre . 122 et Tautre est connectee a la sortie de Tunite 120 par 
rintermediaire d'un organe de multiplication 130 qui multiplie ou non par deux le code de sortie de i'unlte 

40 1 20. Le multiplexeur 125 est muni, lui, de 3 entrees, la premiere est reliee a la sortie de Torgane 110 et les 
deuxieme et troisieme a des organes de decalage 132 et 133. respectivement. Ces organes de decalage 
decalent ou non de 8 bits vers la droite les nombres contenus dans des registres 136 et 137. 

Un registre 140 qui est du type registre a decalage. enregistre ia valeur de I'element binaire de plus 
fort poids du code a la sortie de Tunite 120. Les sorties paralleles de ce registre sont connectees aux 

45 entrees de donnees d'une memoire 145 via un multiplexeur 174 a deux entrees qui permet aussi de mettre 
en relation la sortie del'unite 120 avec la memoire 145. Les composantes de spectre sont rendues 
disponibles a la sortie 150 de la memoire 145, ce qui constitue en fait la isortie du circuit 7. 

Les memoires 105 et 145 sont adressees de la meme maniere que les memoires 61, 63 au moyen d'un 
code d'adresse ADM et ADSP. La memoire 115 qui est une memoire a lecture seule est adressee par le 

50 code ADOS. Les signaux S12, SA2, SB2 et SS2 positionnent les multiplexeurs 107, 124, 125 et 147. Le 
signal M2 applique a Torgane 130 pour autoriser cu non une multiplication par deux des signaux D8 et D8 
force les organes 132 et 133 a diviser par 8 ou non en effectuant un' decalage correspondent a cette 
division. Le code 0PA2 commande les operations a effectuer par I'unlte 120. Le circuit 100 elabore ces 
differents signaux et codes ADM. ADOS, ADSP, SA2, SB2, 0PA2, SS2, D8, D8\ M2 et SI2 en conformite 

55 avec le fonctionnement de ce circuit 7 explique ci-dessus. 



4 - FONCTIONNEMENT DU CIRCUIT D^ANALYSE 7 
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4-A TRANSFORMEE OE FOURIER DISCRETE (TFD) TFD - A1 - Calcul des composantes (phase CF) 

Durant cette phase on effectue le calcul. pour les differentes raies de leurs composantes resile SR(k} et 
imaginaire Sl(k) : 

SR{k) = E X(i).cos(2irk i) (D1) 

i=0 . N 

N-1 - i ■ ■ 

Sl(k) = Z X(i).sin(2irk.^) (D2) 

i=1 . N . 

Les differents echantillons ont ete stockes au prealable dans la nnemoire 105. Le calcul de la 
composante SR(k) comnnence a I'lnstaht tlOO <voir figure 15). A ctiaque coup d'horloge les echantillons X- 
(0), X(1). .... X(N-1) vont venir occuper le registre 112 et simultanement les differentes valours de cos : 1, 
cos277k/N. cos2irk N-1/N, vont venir occuper le registre 113 de sorte que le produit indique par la 
formule (01) puisse s'effectuer. Le resultat du produit est stocke a chaque coup d'horloge dans le registre 
123. Le registre 122 sort d'accumulateur et est mis a zero (INIT) au prealable. Lorsque le calcul de SR(k] 
est termine, le resultat est mis dans le registre 136 (t107). L'instant t107 marque aussi le demanrage du 
calcul de Sl(k) qui s'effectue de la memo maniere que SR(k) mais en commenQant avec rechantillon X(1) et 
les valours de sinus telles que donnees par la formule (D2), Le resultat est mis a TinstantJIII dans .le 
registre 137. 



4-A TFD - A2 - Calcul de la composante S 

Selon une caracteristique de i'invention, la quantite : . 
S(k) = V[SR(k)]^ ~ [SI(k)P (D3) 
est evaluee de la maniere suivante : 
S{k) = MAX[|SR(k)|,|SI(k) + 3/8 MIN[|SR(k)|. |SI(k)|] 
ou MAX[A,B] = Asi A>B 

= B si B>A , 
ou MIN[A.B3 = A si A<B 
= B si B<A 

Cette approximation amene une imprecision de 5 % seulement. Pour cela a I'instant till la valeur SR- 
(k) est transferee du registre 136 au registre 123. Le code 0PA2 est tel que, a la sortie de I'unite 120, on 
obtient sa valeur absolue |SR| qui est transferee au temps t112 dans le registre 122. A cet instant 1112 la 
valeur S!(k) est mise dans le registre 123 de sorte qu'a Tinstant til 3 sa valeur valeur absolue |S!| est mise 
dans le registre 137 et se retrouve a I'instant suivant til 4 dans le registre 123. II est des lors possible de 
determiner la plus grande de ces quantites [|SR| et |SI| appelee f\/IAX et par consequent la plus petite MIN. 
Le resultat de cette comparaison est envoye au circuit sequenceur 100 de sorte que les contenus des 
registres 136 et 137 puissent etre judicieusement envoyes dans les registres 123 et 122. Ainsi, a {'instant 
t115 la valeur MIN qui est contenue dans Tun des registres 136 et 137 est mise dans le registre 123 ; cette 
valeur est divisee par I'unite 120, le resultat est mis a Tinstant t116 dans le registre 122 alors que lea valeur 
MAX est mise dans le registre 122 qui va servir d'accumulateur pour accumuler successivement la valeur 
MAX a rinstant t117 et la valeur MIN/8 a Tinstant t118 obtenue par Tun des decaleurs 132 et 133 ; seion. 
que la valeur MIN est contenue dans le registre 136 et 137 a i'instant tnS la composante spectrale peut 
etre enregistree dans la memoire 145. . 

4-B CALCUL DU SPECTRE A PARTIR DES COEFFICIENTS DU PREDICTEUR 

Les calculs a effectuer sont : . . _ . 
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DR(k) =1+1: a{i).cos2Trk i (D4) 
i=1 ^ 

.5 DI(k) = E a(i).sin2irk ^ (D5) 

D(k) = VDR{kF +"DKkF (D6) 
S(k) = 1/D(k) (D7) 

70 On remarquera que les calculs (D4), (D5), (D6) sont pratiquement identiques a ceux de la TFD el sont 
represenies entre les instants t200 a t228 a la figure 16. La division indiquee par la formule (07) va faire 
intervenir le registre 140. On se place a I'instant t217 ou se termine Texecution de la formule (D5) mais, id, 
le resultat est envoye dans le registre 136 pour se retrouver a Tinstant t218 suivant dans ie registre 123 qui 
contient des lors le diviseur de ta division indiquee a la formule (07) ; en ce meme instant le chiffre "1 " est 

75 mis dans le registre 122 en tant que dividende. Ce chiffre "1" resulte d'une commande 0PA2 appliquee a 
I'instant ant^rieur. Le processus de division demarre alors et le resultat fina{, contenu dans ie registre 140, 
peut alors etre transfers dans la memoire 145 a I'instant t221. 

20 5 - ACCES A L'HORLOGE DE DATATION ET SON REGISTRE 

Un code approprie est applique a I'acces 5 sur ia ligne 13 ; ce code agit par le multiplexeur de sortie 
35 pour mettre en relation I'acces 4 soit avec Thorloge de datation 30 pour une certaine donnee emise sur 
la ligne 10, soit avec Ie registre de datation pour une autre donnee. 
25 On remarquera que le nombre d'acces relativement reduit permet au dispositif de I'invention d'etre 
realise en une seule puce (chip). 

Revendications 

30 ■ . ' 

1. Dispositif d'analyse spectrale d'un signal numerique, caracterise en ce quMi comporte des acces de 
commande pour recevoir une information de methode d'analyse destinee a commander au moins un circuit 
d'analyse. 

2. Dispositif d'analyse spectrale selon la revendication 1, caracterise en ce qu'il comporte un premier circuit 
35 d'analyse pour appliquer une premiere methode sur le signal numerique et pour appliquer au moins une 

deuxieme methode en cooperant avec un deuxieme circuit d'analyse. 

3. Dispositif d'analyse spectrale selon la revendication 2, caracterise en ce que la premiere methode est 
une methode de transformee de Fourier discrete tandis que la deuxieme methode est une methode de 
coefficients de prediction etablis par le deuxieme circuit d'analyse. 

4c 4. Dispositif d'analyse spectrale selon la revendication 3, caracterise en ce que I'information de methode 
peut comporter une sous-information de coefficient pour definir une methode de determination de coeffi- 
cients a etablir par le deuxieme circuit d'analyse. 

5- Dispositif d'analyse spectrale selon la revendication 4, caracterise en ce que la methode de determina- . 
tion de coefficient est Tune et/ou Tautre des methodes suivantes : 
. 45 ' methode de Levinson, 
- methode du gradient 

6. Dispositif d'analyse spectrale selon Tune des revendications 1 a 5, caracterise en ce qu'il comporte. un 
ensemble de registres pour contenir I'information de methode et les parametres qui y .sont associes. 

7. Dispositif d'analyse spectrale selon I'une des. revendications 2 a 6, caracterise en ce qu.e lesdits circuits 
so d'anaiyse ont une structure de processeurs formes de memoires adressables, d'unites de calculs et 

d'organes de sequencement. 

8. Dispositif d'analyse spectrale selon I'une des revendications 1 a 7, caracterise en ce qu'une correlation 
du signal numerique est effectuee et en ce qu'une sortie est prevue pour fournir le coefficient R(D) de la 
correlation. 

55 9. Dispositif d'analyse selon I'une des revendications a 8. caracterise en ce que I'analyse de spectre 
indique une operation du type : 

D(k) = v'DR(k)^ ~ DRl^p 

et en ce que le circuit d'analyse compcrtt?. des premiers moyens pour determiner la valeur MAX de |DR(k)| 
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et de |DI(k)| et la valeur MIN de |DR(k)| et de |DI(k)| des deuxierhes moyens pour effectuer I'operation : 
D(k) - MAX(|DR|.1DI|) + 3/8 MIN{|DR|,|DID. 

10. Dispositif d'analyse spectrafe selon Tune des revendications 1 a 9, caracterise en ce qu'il est realise en 
une seiile puce. 

5 11. Dispositif d'analyse spectrale selon Tune des revendications 1 a 10, caracterise en ce qu'il comporte 
une horioge de datation accessibre a Tutilisateur pour dater certaines donnees. 

12. Dispositif d'analyse spectrale selon la revendication 11, caracterise en ce qu'il comporte un registre 
pour stocker la date fournie par I'horioge de datation. \ - 

TO ■ . • 
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